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For future surveys,0^^ spectroscopic follow-upiifor all supernovae will be extremely diffi- 
cult. However, one can use light curve fitters, eIS to obtain the probability that an object 
is a Type la. One may consider applying a probability cut to the data, but we show 
that the resulting non-la (nia) contamination can lead to biases in the estimation of 
cosmological parameters. A different method, which allows the use of the full dataset 
and results in unbiased cosmological rparameter estimation, is Bayesian Estimation Ap- 
plied to Multiple Species (BEAMS). B BEAMS is a Bayesian approach to the problem 
which includes the uncertainty in the types in the evaluation of the posterior. Here we 
outline the theory of BEAMS and demonstrate its effectiveness using both simulated 
datasets and SDSS-II data. We also show that it is possible to use BEAMS if the data 
are correlated, by introducing a numerical marginalisation over the types of the objects. 
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An Introduction to BEAMS 

Consider the following thought experiment: there is a room full of people and we 
desire to know the average height of women in the room. However, the only informa- 
tion we have is the height of each person and the length of their hair. Now, it is easy 
to imagine we could write down the probability that a person is a woman based 
on their hair length. How, using those probabilities can we estimate the average 
height of women? We could try applying a probability cut, for example assuming 
that anyone with a probability higher than 0.9 is a woman. This approach has some 
obvious problems: not only would we cut out much of the data (women with shorter 
hair) but also men with longer hair will be included and will bias the estimate of 
the average height. 
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This thought experiment is directly analogous to supernovae cosmology: we want 
to estimate cosmological parameters, such as Qrn and Hq, from supernova data. 
We can obtain probabilities from the light curves of the supernovae, but can only 
be sure of the type with spectroscopic follow-up. Only the la's give cosmological 
information, like the women in the sample, but we cannot ignore possible contami- 
nation from the nia's, as we could not ignore the long-haired men in the example. 
Bayesian Estimation Applied to Multiple Species, or BEAMS,EI utilises all available 
data, correctly handles the contamination from the tall men and produces unbi- 
ased results. BEAMS includes the uncertainty in the types by marginalising over 
all possible combinations of types in the posterior: 

Pi0\D) = Y,PiO,r\D), (1) 

r 

where is the set of cosmological parameters, D is the data and r is a length 
vector of types (where N is the number of datapoints). The problem with this 
analytic marginalisation is that it is (assuming two types, la and nia), an order 
2^ calculation, which is computationally unfeasible. Fortunately, if the data are 
uncorrelated, Ref. || showed that it is possible to simplify the calculation to order 
2 A". The resulting calculation is essentially a kind of mixture model (see Ref. |7| for 
an excellent video introduction). 

Figure |l| compares BEAMS with a probability cut and with the spectroscopic data 
only, when using both a controlled, simulated dataset and the SDSS-II dataset.@l3 
The SDSS-II dataset is small, so no bias can be seen when a probability cut is used. 
However, as future surveys such as the LSST increase the number of supernovae 
dramatically, we expect a bias to become visible, similar to the constraints shown 
in panel (a). 
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(a) (b) 

Fig. 1. Qrn — contours produced by analysing (a) simulated data and (b) SDSS-II data, using: 
only the spectroscopically confirmed objects (dashed line), a sample with a probability cut of 0.9 
(solid line) and BEAMS (filled contours). The filled square, empty square and the cross-^iepresent 
the fiducial cosmology, the BEAMS best fit and the spectroscopic best fit respectively. t3 
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Photometric supernovae data can have many sources of correlated systematic 
uncertainties, such as filter errors, peculiar velocities and redshift-stretch correla- 
tions.ll2l If the uncorrelated form of BEAMS is used to analyse this data, biases can 
result. We developed a methodic of using BEAMS for correlated data by marginal- 
ising over the types numerically instead of analytically. Figure || shows an example 
covariance matrix used (a) with the corresponding Qrn — ^A contours (b), for a mock 
dataset of 200 strongly correlated supernovae. It is clear that without accounting 
for the correlated systematic uncertainties, it is possible to obtain biased contours 
which severely underestimate the errors on the parameters. 
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Fig. 2. (a) An example covariance matrix,0 where only the la's are correlated (and the diagonal 
has been removed for clarity), (b) flm — contours produced when analysing the mock using the 
uncorrelated form of BEAMS (grey, dotted) and the correlated form with numerical marginalisation 
(red, filled). 
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